Kinetic theory of fluctuations in conducting systems 
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We propose an effective field theory describing the time dependent fluctuations of electrons in 
conducting systems, generalizing the well known kinetic theory of fluctuations. We apply then the 
theory to analyze the effects of strong electron-electron and electron-phonon scattering on the statis- 
tics of current fluctuations. We find that if the electron-electron scattering length is much shorter 
than the transport mean free path the higher cumulants of current are parametrically enhanced. 



X3 

o 

> 

m 
o 

o 



I 

o 
o 



X 



Fluctuations of electric current in time have been the 
subject of fundamental and applied research since the 
time of Schottky. They have been a centerpiece in the 
study of systems ranging from vacuum tubes to quan- 
tum Hall bars (for a review see [J). More recently noise 
has been recognized as an important diagnostic probe 
for the nature of disorder and interactions. Tradition- 
ally the number of electrons taking part in transport 
in macroscopic systems is large; for that reason current 
fluctuations obey nearly Gaussian statistics. Under such 
circumstances it is the power spectrum which is the ob- 
ject of primary interest when addressing the fluctuations. 
However, as the sample size shrinks down from a vacuum 
tube to a nanojunction, and the currents involved become 
correspondingly smaller, the central limit theorem does 
not apply any more, and information concerning higher 
current cumulants may be needed. Various techniques 
have been applied to address this issue 0, 0, fllM E| • 

There is an interesting and relevant class of systems 
out of equilibrium in which the effect of quantum inter- 
ference on transport quantities is negligible. Transport 
through such systems is often described by classical ki- 
netic equations (e.g. the Boltzmann equation). The in- 
tensity of the noise spectrum can then found employing 
a related approach, known as the kinetic theory of fluc- 
tuations (KTF) 0,1,13. The appeal of this approach 
is its intuitive transparency, yet traditionally it was con- 
structed to describe pair correlation functions, hence was 
limited to finding the second current cumulant. In a re- 
cent insightful work Nagaev T^l has proposed a gener- 
alization of the KTF, expressing higher order cumulants 
in terms of pair correlators. Using this idea (which was 
termed the "cascade approach"), Nagaev calculated the 
third and the fourth cumulants of the current. While this 
approach was shown to successfully address a number of 
problems, its (quantum-mechanical) microscopic founda- 
tions, as well as its formulation in terms of an effective 
field theory (which would allow for further generaliza- 
tion, such as including interactions or calculating the full 
counting statistics (FCS)) were not fully understood. An 
important step in the latter direction was an introduc- 
tion of a stochastic path integral by Pilgram, Jordan, 
Sukhorukov, and Biittiker [lll |. 

In the present paper we propose a phenomenological ef- 



fective field theory describing the time dependent fluctu- 
ations in phase space of electrons in conducting systems. 
This is a field-theoretical implementation of the Nagaev's 
regression scheme. This theory allows for various gener- 
alizations, including the effects of electron-phonon and 
electron-electron interactions. In a number of instances, 
this theory is equivalent (up to small quantum correc- 
tions) [T^ IT^Ito the non-linear a-model formulation of 
the problem [g. In particular, for non-interacting elec- 
trons it reproduces the FCS ^14l| of current fluctuations 
obtained earlier within quantum approaches. We apply 
the theory to systems with a strong electron-electron 
scattering ("hot electron regime"). We find that (in 
the absence of umklapp processes) , the current statistics 
changes dramatically depending on the relation between 
the electron-electron scattering length Zc-c and the elas- 
tic mean free path. Specifically, in the regime ^o-o <C hr, 
when the electrons form a weakly disordered liquid, the 
current cumulants are parametrically enhanced due to 
the electron-electron scattering. We find the scaling law 
governing the dependence of the FCS of hot electrons on 
the ratio Finally, we show that under strong 

electron-phonon interaction the statistics is Gaussian. 

The effective action. The KTF can be cast as a Boltz- 
mann equation with a Langevin-type noise 0, 



£ + / /(p,x,i) = W(p,x,i). (1) 



Here we have introduced the Liouville operator 

ot m ox op 

while /(p, X, t)(ix(ip is the number of electrons in a 
unit cell centered at (p, x) in a phase space at time 
t. With p(x, i) = e /(p, x, t) being the local den- 
sity of electrons, the self-consistent electric field satisfies 
divE = Anp. 1 The operator / = /dis + lo-ph + ^o-o 
is the collision integral Q with the respective con- 
tributions coming from disorder, electron-phonon and 
electron-electron scattering. The Langevin random 
source SJ{p, x, t) in ^ represents the fluctuations of the 
net flux of particle in phase space at (p,x). These fluc- 
tuations are the consequence of the scattering processes 
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described by the respective collision integrals, and should 
be treated consistently along with the collision integrals 
to account for the conserved quantities. For the physics 
considered here (and disregarding quantum corrections) 
the interplay among the various scattering channels can 
be neglected. Within the same accuracy the Langevin 
flux is the sum of three uncorrelated random processes 
5 J = SJdis + SJee + SJo-pii. Its correlation function can 
be found employing semiclassical considerations . The 
random flux correlator can be written schematically 
as {SJ{p, X, t)SJip',x',t) = (5(x - x')Sit - t')H[p, p', x], 
where H[p,p',x] = i?dis[p, p', x] + i/eo[p,p',x] + 
iJe-ph[p, p', x]. Explicit expressions for the collision in- 
tegrals and the pair correlation functions of the random 
fluxes are available 

We consider a two-terminal system subject to an exter- 
nal voltage bias V; the ambient temperature at the reser- 
voirs is T. Transport through this system is a stochastic 
process. The charge transmitted over a time interval t 
is Q{t) — dt'I{t') where I{t) is the electric current. 
The number of particles n — Q/e is a, fluctuating quan- 
tity described by probability distribution function P{n). 
The Fourier transform of the latter k(A) = ^„ e''^"P(n) 
is the characteristic function which provides a complete 
statistical description of the problem. Constructing a 



field theory it is natural to look for a characteristic func- 
tion in a functional integral representation Q 

k[\] = J Vfip, x, t)VSf{p, X, t) exp(*5[/, Sf; A]). (3) 

Here f,Sf are independent bosonic fields and the mea- 
sure of integration is flat. The field /(p, x, t) is the phase 
space density operator within the KTF; A(p,x, t) is an 
auxiliary 'counting field'. The field Sf{p,x,t) accounts 
for fluctuations of this 'counting field'. The effective ac- 
tion consists of two parts 

^[A] = 5Dy„ + 5Noisc[A]. (4) 

The propagation of fluctuations in phase space is de- 
scribed by the dynamical part of the effective action: 

oo 

iSuyn = i d,yi dt-^Sf{p,x,t)(c + ijf{p,:x:,t). 



The second part of the action accounts for the generation 
of fluctuations: 
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Sfip',^,t) 



(6) 



The effective action, Eqs. Q and ©, can be obtained by exponentiating the equation of motion ^ with a help of 
the Lagrange multiplier Sf{p,x,t) and averaging the resulting e*'^ over the Langevin fluxes. The counting fleld, as a 
function of p, x, and t should be chosen accordingly to the problem one is interested in. It is illuminating to note that 
once the fluctuations of counting field are neglected in S'Noisc, the effective action, (0} is equivalent to the standard 
KTF (i.e., it yields correctly the current-current fluctuations but not higher cumulants). 

So far our treatment has been general. We now assume that the electron gas is degenerate and consider a quasi 
one-dimensional diffusive conductor, where the momentum relaxation length is much shorter than the system length 
{Iti- L). The problem then simplifies considerably. Additional important parameters are the diffusion coefficient, 
D, the single particle density of states ly at the Fermi level and the cross section S. Together they define the Ohmic 
conductance G = e^vDS/L. Within the diffusive approximation, © takes the form 
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Sf{e,x,t) 



dt 
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de 



Ic-ph + Ic-c >fie,x,t) 



(7) 



where e = /2m is the kinetic energy. Eq.(|7I) describes propagation of fiuctuations in diffusive system in the presence 
of a long range Coulomb interaction. In the absence of the Coulomb interaction, one should set E = 0. 

The noise-generating part of the action will also have three contributions, corresponding to elastic (impurity), 
electron-phonon and electron-electron scattering. The noise generated through elastic scattering is given by 



iS\ 



Noise 



dxdedtD 



X{t) + eE 



df{e,x,t) dSf{e,x,t) 



de 



dx 



f{e,x,t){l~ f{e,x,t)) 



(8) 



One can check that the self-consistent field E does not affect the low frequency current statistics. Indeed, 
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a static electric field E = —dx<p{x) can be eliminated 
from Eqs. Q, l|HJl by redefining the energy, e + e4i{x) e. 
This implies that the counting statistics is the same in the 
interacting and non-interacting problems for frequencies 
much smaller than the inverse Thouless time, ui <C D/L^. 
We will consider this regime below and thus set E = 0. 

We are interested in the statistics of the number of 
particles traversing the system during a certain time win- 
dow, irrespective of how this flow is distributed over that 
interval. We thus select the counting field to be 



A, <t <t . 
0, otherwise. 



(9) 



Here A is the argument of the characteristic function 
«(A). 

To analyze the counting statistics we use a saddle point 
approach, which is justified by a large value of the dimen- 
sionless conductance g = G/{e'^/h) ^ 1. Note that this 
is the same condition which also allows us to neglect the 
quantum corrections to the current statistics due to the 
weak localization and the interaction (Altshuler-Aronov) 
effects. The boundary conditions for the distribution 
function are /(e, 0) = /^(e), /(e, L) = /^(e) where /^(^) 
are the equilibrium distribution functions of the electrons 
on the left (right) reservoir. The auxiliary field Sf van- 
ishes at the boundaries. 

Current statistics and inelastic processes. When the 
inelastic processes can be neglected (which requires that 
the shortest inelastic length is longer than the sample 
size, /in ^ L), the saddle-point equations following from 
O, permit an analytical solution 0, 0, |23|. In 
particular, in the shot-noise limit {eV ^ T) the result- 
ing action reads (including a source term for the average 
current, iS^ = iGVXt) 



iS*^'" = fGFarcsinh^(Ve*^-^ - 1) 



(10) 



The found characteristic function, k(A) = expjiS''*"'} 
ea. (|10|l . agrees with the known result for the counting 
statistics obtained earlier |^ |^ Q| within the quantum 
approaches. This demonstrates the equivalence up 
to 1/g corrections, between our effective theory and the 
cr- model formalism jQ]. Indeed, in this case the above 
effective action can also be obtained from a particular 
parametrization of non- linear ct- model |13] . 

In the opposite limit, li^ <^ L, the collision integral 
plays a major role. There is only a small number of 
hydrodynamic modes (determined by the conservation 
laws) that propagate for a distances longer that hn- For 
electron-electron collisions these are the fluctuations of 
the electron density and the energy (or, in other words, 
of the chemical potential and the temperature). By the 
same token, one has to keep only these two modes in the 
noise part of the action, S'noIsc, Eq- Q- In the pres- 
ence of a strong electron-phonon scattering, /o-ph <C L, 
the collision integral le-ph conserves only the number 
of electrons. The role of the electron-phonon scattering 



thus amounts to projecting the action onto the mode of 
the density (chemical potential) fluctuations. 

We consider now the case of a strong electron-electron 
scattering, Icc <Si L, in detail. After projecting the fluc- 
tuating fields /(p, X, t) and Sf{p, x, t) in lO, (01 onto the 
relevant subspace, we are left with the fields fj, and T that 
parametrize the fluctuations of a chemical potential and 
the temperature, respectively, and the conjugate flelds p, 
and T. In terms of these fields the action (O, ^ takes 
the form 



^S'Dyn = iv I dxdt 
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■D*{T)dTdT'^ + Ddti{dfl + fidT) 



(11) 



i-SNoisc = —V I dx I dtT 



-oo 
.^2 



D{\ + dil+iidTf 



+D*{T)—T'^{dff 



(12) 



Here we introduced a thermal diffusion coefficient 
D*{T) = vUr,-' + r;-i(T))-V3, where 



1 



dpdp nn iJcc[p,P 



T T 



(13) 

with n, n' being the direction of the momenta p, p' on the 
Fermi surface. Physically, is the relaxation rate of 

the energy current due to the electron-electron collision 
processes. As we are going to show, there are two strik- 
ingly different cases. For Ztr <SC leo, which is the commonly 
considered hot-electron regime, the thermal diffusion and 
particle diffusion constants coincide, D* ~ D. Our action 
((TT|l . ((T^ then agrees with that obtained in 
framework of a stochastic path integral formalism. In this 
regime, the electron-electron scattering does not affect 
the FCS qualitatively but only modifies the numerical 
coefficients in the expressions for current cumulants 
For example, for el/ > T one finds ^2 = {V^/'i)e'^GV , 
S3 = (8/7r2 - 9/16)e^GV, etc., instead of ^2 = ^e^GV, 
S3 = -j^e^GV, ... for non-interacting electrons. 

The situation is qualitatively different in the regime of 
very strong electron-electron scattering (or, alternatively, 
of a very clean system), l^e <C ^tr, when the electrons form 
a liquid flowing through a weak disorder potential. In 
the absence of umklapp processes, the electron-electron 
scattering conserves the total momentum, so that the 
density diffusion is much faster than the thermal diffu- 
sion, = D/D*{T) 3> 1. The saddle-point equations 
following from (|ll|l . (|12|l and determining the FCS of 
hot electrons do not permit an analytic solution. To ob- 
tain a qualitative understanding of the problem we first 
consider a simplified model, approximating D*{T) by a 
constant D* . In this case, the action (|ll|l . (|12|l allows 
us to determine the scaling dependence of the FCS on a. 
By performing a rescaling of fields fl — > fl/a^ fi ^ fi/a 
in l|ll|l . H12|l we find that k(A) has the following scaling 
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form, 

k(A) = cxpi-a-^GTtJ'iaX, aeV/T)}, (14) 

where the function !F{x,y) governs the FCS at a = 1, 
k(A) = exp{-GTt J^{X, eV/T)}. This impHes that in 
the equihbrium hmit, T ^ aeV , the cumulants Sn with 
n > 3 in the electronic hquid are strongly enhanced com- 
pared to their values in the conventional hot-electron 
regime (a = 1). Specifically, the enhancement factor 
is a^™ for the cumulants S2m+i and 52rn+2- In the op- 
posite limit (r — > at finite V), the enhancement starts 
already from 52, with 5„ enhanced by the factor a"~^. 
Physically, the enhancement of higher moments of cur- 
rent fluctuations is due to the coupling between the mode 
of electric conductivity and the slow mode of the thermal 
conductivity. 

So far we have used a model of a temperature in- 
dependent thermal diffusion coefficient. On the other 
hand, in the Fermi-hquid theory, D*{T) ~ epVp/T'^. 
We argue now that in a realistic situation the above 
strong enhancement of the cumulants in the ^oo ^ hr 
regime remains parametrically valid, with replaced by 
D* / D(Tcfi), where Toff is the characteristic temperature 
in the system. In particular, in the equilibrium situation 
(V^ = 0), Teff is simply given by the temperatute T of the 
reservoirs, ~ Tti/Tee{T). In the non-equilibrium limit 
(T = 0, y finite), the effective temperature is found to in- 
crease exponentially with V TS] , yielding asymptotically 
~ exp{{3/4TT^)n,/T,,{V)} 

Finally, we consider the case of a strong electron- 
phonon scattering. The only soft modes allowed in this 
case (^c-ph ^ L) arc fluctuations of chemical potential. 
The saddle point equations d^p, — , d^n — along 
with the boundary conditions lead to the trivial solution 

fi{x) = m(0) + ix/L)eV, fi{x) = 0. (15) 
The resulting action 

■^.c-ph ^ „fG[TA2 + iVX] (16) 



represents Gaussian fluctuations, so that higher cumu- 
lants are suppressed by the electron-phonon scattering. 
This is consistent with earlier predictions based on qual- 
itative arguments T| and calculations of the third cumu- 
lant 6]. However, to the best of our knowledge, it is the 
first analytic derivation of this general result. 

To summarize, we have proposed an effective theory 
describing the creation and evolution of fluctuations in 
phase space for degenerate electron gases and liquids. We 
have employed this theory to analyze statistical proper- 
ties of noise in disordered conductors. We have repro- 
duced known results for free fermion and conflrmed the 
prediction of Gaussian statistics for the case of strong 
electron-phonon scattering. We have analyzed the effect 
of electron-electron scattering on the statistics of current 
fluctuations and established a scaling relation governing 
the dependence of the FCS on the ratio hr/lee, show- 
ing that the statistics differes parametrically between the 
^tr lee and hi ^ lee cascs. Specifically, in the latter, 
when the electrons form a liquid fiowing through weak 
disorder, the current cumulants are parametrically en- 
hanced by electron-electron scattering. This is related 
to the coupling of the density mode to the slow mode of 
temperature fluctuations. The regime Itr ^ lee is thus 
most favorable for experimental investigation of higher 
cumulants of current. 
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Note added. As this work was nea ring co mpletion we 
have learned of several related works [13, uB, l20| that 
appeared in the cond-mat preprint archive. 
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